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Abstract. A Green's function method is developed for solving strongly-coupled 
gravity and matter in the semiclassical limit. In the strong-coupling limit, one 
assumes that Newton's constant approaches infinity, G ^ oo. As a result, one 
^ may neglect second order spatial gradients, and each spatial point evolves like an 

ly-^ ■ homogeneous universe. After constructing the Green's function solution to the 

O^l . Hamiltonian constraint, the momentum constraint is solved using functional methods 



o 

(N 



in conjunction with the superposition principle for Hamilton-Jacobi theory. Exact and 
' approximate solutions are given for a dust field or a scalar field interacting with gravity. 

00 ■ 

O . Strong-coupling expansions have proven to been useful in elucidating gauge theories 

^ ■ of modern particle physics (see, e.g., Creutz 1983 and Polyakov 1987). A general class 
bX)' of Hamilton-Jacobi solutions for strongly-coupled gravity and matter will be presented 
here. It is hoped that these solutions will help clarify quantum features of gravity 
^ ■ that are necessary in formulating the inflationary scenario (Hartle & Hawking 1983 and 
Hartle 1997). 

In the context of classical gravity, the strong-coupling expansion is very similar 
although not exactly the same as a long-wavelength expansion. In their analysis of 
singularities in the early Universe, Lifshitz and Khalatnikov (1964) expanded the 3- 
metric in a Taylor series in what was basically a strong-coupling expansion (see also 
Belinski et al 1970 and Landau and Lifshitz 1975). However, their program was 
incomplete in that they did not give an explicit solution of the momentum constraint 
equation which is just the Einstein equation. Pilati (1982), Teitelboim (1982) and 
Henneaux et al (1982) took some initial steps in formulating a quantum treatment of 
strongly-coupled gravity although many difficulties still remain (see also Husain 1988). 
They pointed out that this system is much simpler than the case for full gravity: the 
Hamiltonian was similar to that of an homogeneous minisuperspace model. Solutions 
were given which were a product over the points on a lattice. Once again, they did not 
attempt to solve the momentum constraint. In the present work, I hope to improve upon 
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this situation by showing how to simultaneously solve both the energy and momentum 
constraints in the semiclassical limit. Special cases of the semiclassical limit were given 
before (Salopek 1991), but now general solutions will be given. 

A long-wavelength expansion has been considered by numerous authors, including 
Tomita (1975), Salopek and Stewart (1992), Deruelle and Langlois (1995) and Veneziano 
(1997). However, by using Hamilton- Jacobi (HJ) methods, one can elegantly solve the 
momentum constraint. In fact. Parry et al (1994) used a spatial gradient expansion 
which is manifestly invariant under spatial coordinate transformations to give a HJ 
solution for gravity and matter. They described a systematic method of obtaining 
terms of arbitrarily high order. Soda et al (1995), Chiba (1995) and Darian (1997) 
have extended this method to encompass Brans-Dicke gravity, n-dimensional gravity, 
and gravity interacting with electromagnetism, respectively. 

In a previous paper, Salopek (1991) had shown that one could construct solutions 
of the energy constraint and the momentum constraint for strongly-coupled gravity by 
assuming that the generating functional was of the simple form: 

5 = -2 1 d'x^'/' H[(l)ix),xix)], (1) 

where 0(x) and xi^) denote the matter fields. If the Hubble function H were a constant, 
then this functional would correspond to the volume of any given 3-geometry. For 
many astrophysical applications, the ansatz of eq.(|l|) is sufficient. Such a formalism 
provides an elegant description of the long-wavelength evolution of fluctuations arising 
from inflation. For example, it was used to construct inflation models that yielded 
non-Gaussian fluctuations (Salopek 1992a). Such models are still of observational and 
theoretical interest (Moscardini et al 1993). However, for future extensions of the HJ 
approach, it is of interest to construct a more general class of solutions. Some steps 
were taken in that direction by Salopek (1991) who presented a complete solution of 
the energy constraint of the long-wavelength problem. In addition, by a suitable choice 
of the spatial coordinates, one could solve the momentum constraint. In the present 
work, powerful geometric methods will be utilized to give a more elegant solution of the 
momentum constraint. 

Over the past few years, an ever increasing number of tools have been developed in 
order to solve the functional Hamilton- Jacobi equation. They include: 

(1) The Spatial Gradient Expansion. For centuries, Taylor series solutions have been 
constructed for linear equations. Series expansion techniques can also be applied to 
nonlinear differential equations including the HJ equation for general relativity (Parry et 
al 1994, Salopek 1997). In this way, it is possible to decompose semiclassical superspace 
into a sum of minisuperspaces. 

(2) The Superposition Principle. Since the equations for quantum mechanics are linear, 
the superposition principle plays an important role. By employing the stationary phase 
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approximation, one can also enunciate a Superposition Principle for Hamilton- Jacobi 
theory (see, e.g., Salopek 1997). Although it is not hnear, it proves very effective 
in constructing complicated solutions from more elementary ones (Landau & Lifshitz 
1960). It will also prove useful in solving inhomogeneous gravitational problems such 
as the one presented in this paper. 

The Lagrangian and Hamiltonian formulations for strongly-coupled gravity and 
matter are given in section 1. The semi classical solution method may be simply 
explained using some rudimentary ideas from quantum mechanics. A complete solution 
of the energy constraint for a dust field interacting with gravity is described in section 2. 
This solution is related to the well-known Kasner metric. The momentum constraint is 
solved using a superposition over the complete set of solutions. It is also shown how to 
construct the set of constant junctionals. Explicit solutions are constructed for a single 
dust field in section 3 and for two dust fields in section 4. In sections 5, 6 and 7, the 
entire analysis is repeated for a scalar field interacting with gravity. Conclusions are 
given in section 8. 



1. Strongly-Coupled Equations for Gravity and Matter 

The action for Einstein gravity interacting with a dust field, x, and a scalar field, 0, 
may be written as follows: 

jd'x^ i^l^^)R-l-^g^^4,^^4,^,-nV{ct^)-''-^^^^ ,(2a) 

where 

K = SttG = Sn/ml {2b) 

is the gravitational coupling constant, n = n{t, x) is a Lagrange multiplier which ensures 
that the square of the four-velocity, U^, is minus one: 

U^^-g^''x,u/i^. (3) 

The above form can be obtained from the usual one by appropriate scalings of the matter 
fields and their coupling constants. Factors of k appear in key places of eq.(2a) in order 
to obtain the desired energy constraint: 

-1/2 
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1/2 1/2 

K 7rVl + xl«X|a/«^ - V ^ + ^<P%a + « l'^' ^(0) = . (4a) 



2k 2k 

No factors of k appear in the momentum constraint: 



Hiix) = -2 (^ikTr"') . + Tr'^SfcM + ^% + 7r^X,i = . (4&) 



In a classical Hamilton-Jacobi formulation of general relativity, one defines the 
generating functional, 

S = S[yab{x), (l){x),x{x)], (5a) 

by assigning a real number to each field configuration [^(a;), x(a:)] on a space-like 

hypcrsurface with 3-geometry given by 7a6(ic). In a semiclassical context, one allows 
for the possibility that S may be complex. The Hamilton-Jacobi equations are obtained 
from the constraints (4a-b) by replacing the momenta with functional derivatives of S: 

7r^{x) = -^, Ax) = -^ n^\x) = -^ {5b) 
dx{x) d(f){x) S-fab{x) 

In the limit of large gravitational coupling, k — > oo, the generating functional S^^^ 
for the strongly-coupled system evolves according to the following equations: 

n^^\x)/K^ r'l' (27ac76.-7a.7c.)^^+ ^ 

"7a6 "7c(i ''X 



+ r;^!f^f!lV^y/Vw=o. (6a) 



The energy constraint is ultra-local in the sense that different spatial points are not 
coupled to each other. As a result, the Poisson bracket of 1-L^^\x) with 7i''^\y) vanishes, 
and consistency of the Hamiltonian constraint at different spatial points is assured. 
Consistency or 'integrability' of the Hamiltonian constraint for full general relativity is 
more complicated. It is related to the freedom in choosing an arbitrary time foliation: 
see Parry et al (1994) as well as Salopek (1995). However, the momentum constraint 
does couple the spatial points together because spatial derivatives appear. It is this 
cross-coupling of different spatial points that makes the strongly-coupled system non- 
trivial. 

(In order to simplify the notation, S will be used to denote S^^^ for the remainder 
of this paper.) 



1.1. Analogy from Elementary Quantum Mechanics 

The general solution for the strongly-coupled system may appear quite complicated so 
it is instructive to consider a well-known example from elementary quantum mechanics 
which illustrates the essential features. 
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1.2. Rotationally Symmetric Solutions to the Schrodinger Equation for a Free Particle 

The two-dimensional Schrodinger equation for a free particle with zero angular 
momentum is given by: 

.9^ 1_ &^ 

dt 2m \ dx^ dy"^ j 

In solving these two equations, one would ordinarily invoke polar coordinates, (r, 9), and 
then write the solution as ilj = ijj(r). Instead, 1 will utilize a circuitous method which is 
of interest because the same technique may be generalized to the case of strongly-coupled 
gravity. 

One notes immediately that the first equation (7a) can be solved generally without 
any reference to the second by using a superposition of plane waves: 

^ = j (fk f{k) e*(^-^^-'"*) w={kl + kl)/{2m), (8) 

— * — * 

where f{k) is an arbitrary function of the wave-vector k. Please note that a plane wave 
^t{k-x-u}t) jg ^ solution of the symmetry condition, eq.(7&). However, by suitably 
restricting the function / so that it is function only of the magnitude of A; = \Jki-\- /cl, 

f{k) ^ f{k) (9) 

one can indeed satisfy the symmetry condition. 

In a scmiclassical gravitational context, the energy constraint (6a) will be analogous 
to the Schrodinger equation (7a) and the momentum constraint (66) will be analogous to 
the symmetry condition (7&). In general relativity, it is not yet known how to explicitly 
solve the momentum constraint such that one obtains a set of reduced variables {physical 
variables) analogous to solving the constraint and then deducing that is solely a 
function of r. Instead, one solves the strongly-coupled system by first obtaining a general 
class of solutions for the energy constraint. By considering a suitable superposition over 
these solutions, one then constructs a solution to the momentum constraint. 

Shortly after Dirac (1958) gave general relativity its Hamiltonian form, Higgs 
(1958) pointed out that the momentum constraint of general relativity implied that 
the wavefunctional was invariant under spatial coordinate transformations. Although, 
the Hamilton-Jacobi equation was written explicitly by Peres (1962), only quite recently 
(Salopek and Stewart 1992, Parry et al 1994) has the spatial coordinate invariance been 
exploited to construct explicit solutions of the full HJ equation. They utilized a spatial 
gradient expansion. 
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2. Strongly-Coupled Solutions: Dust Interacting with Gravity 

In the strongly-coupled limit, the energy constraint for gravity and a dust field is: 

O'Jab O'Jcd OX 

2.1. Green's Function Solution for the Strongly- Coupled HJ Equation 
A complete solution, Q, to the above equation is: 



(10) 



G [lab{x) ,X{X)\ lab (x) , X(0) (x)] 



27^^^ lii) coshUl ^z) -7 



1 



{x{x) - X(o){x] 



1/2 



1/2 



where 



^^Tr{ln([/.][/.(,5])ln([/.][/.(,;])} 



and [h] and [h^^^] are matrices with components given by 



[h]ab = l ^^^lab, 



;ila) 



116) 



:iic) 



This solution will be referred to as the Green's function. It depends on 7 parameter 
fields, X(o) (x) and 7^°'' (x) , which may be interpreted as the initial fields for the dust and 
3-metric. Although the variable z had been defined in an earlier paper (Salopek 1991), 
the present form for the generating functional Q is new. It is superior to the previous 
formulation in that it is symmetric upon interchange of the metric jabix) with the initial 
metric 7a°'*(a;): 

Gbab{x),x{x) \ labix)^X(0){x)] = Qllab i^) ^ X{x) \ lab{x) , X{0){x)] . (12) 



One may verify by differentiation that eq.(lla) is indeed a solution to eq.(|10]) but in its 
present form it does not satisfy the momentum constraint. 



n^\x) = -2{r. 



6S \ 6S 6S 

'i^T — +7 — lki,i + ^ X,i = , 

oikjjj hki ox 



(13) 



because the parameter fields are spatially dependent and hence they carry momentum. 



2.1.1. Too Many Parameter Fields There is an additional problem with the Green's 
function solution eq.(lla) as it stands: there is one too many parameter fields. This 
subtle fact is not immediately apparent but it becomes painfully clear after one explicitly 
tries to implement the superposition principle in Section 2.2. One can resolve this 
problem by arbitrarily setting one of the parameter fields to zero, say X{o){.x) = 0. A 
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brief explanation of this point now follows although it is best to skip over the next two 
paragraphs upon first reading of this paper. 

For each spatial point, there were originally seven degrees of freedom associated 
with x{^) ciiicl the symmetric 3- metric jabix)- Because of the symmetry between the 
initial fields, [j^b {x), X{o){x)] and the original fields, there is an equal number of degrees 
of freedom associated with the former. I will now demonstrate that the initial fields are 
not really independent. First note that there is a relationship between 



because the Hamiltonian constraint is also valid when expressed in terms of the initial 
fields: 

= n^^\x)/K = 7(o;/' (27i^J? - 7l?7l?) vr^o') ^i^) + ^fo) ■ (15) 

If one were to apply the superposition principle in Section 2.2 by minimizing Q + Q with 
respect to j'^Jix) and X{o){x), one would find that 

As a result, the momenta 

'<°)<^' = ^^' 

are functions of 'j^^ix) and X{o){x)- Hence eq.(|T5]) implicitly defines a relationship 
amongst the initial fields [j^b {x),X(o){x)]'- they are not independent variables. 

In order to avoid the above problem I will now arbitrarily fix one of the parameter 
fields. In the present work, I will choose X{o){^) = 0, while 'j'^b{x) will remain a 
free variable. The subsequent analysis is indeed consistent, although other choices are 
presumably possible. 



2.2. Superposition Principle for HJ Theory and Solution of the Momentum 
Constraint 

Assuming X{o){x) = in the general solution to the Hamiltonian constraint, eq.(lla-c), 
one can construct a solution S to the momentum constraint through a superposition 
principle: 

S[^abix),xix)] = gbabix),xix)\ (^) , X(o) (x) = 0] + Q^l?] , (18a) 
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where 'y^^ix) has been chosen to minimize G + Q, 

hah hab 

Here Q is an arbitrary "gauge-invariant" functional of the initial 3-metric, 7^^'': 



(0). 



(0) \ (0) 



Construction of solutions to the strongly-coupled problem through the superposition 
principle eq.(18a-c) will be known as the Semiclassical Green's Function Method. 

An elementary discussion of the Superposition Principle for HJ theory was given by 
Salopek (1997). As will be demonstrated in section 3.3, the functional Q[7a°''] niay be 
interpreted as the initial state for S when x{^) — 0" 

S[^ab{x),x{x)^0]^Qbab]. (18c/) 

Since the Poisson bracket of Ti,l'^\x) with Ti^^^x) vanishes (weakly), the momentum 
constraint is preserved upon evolution: if S satisfies T-lf'\x) = for x(a;) = 0, it is 
guaranteed to solve it for arbitrary x(a:). 

2.3. Classical Evolution 

Classical evolution is obtained from the minimization prescription eq.(186) by solving 
for ^ab{x) in terms of 'y^ix) and xi^)' 

- (1 - ^x7ro;%))' - 1 x%] ^,^^2 , (i9«) 

. = /ftanh- {^x (7(0} ^S^'^ / (l " ^ Xlo'^'-o)^ , 



[h] = exp 



2z 



TT 



(0)_ 



1^(0) ^ob ) 



(0)\V2 



;i9c) 



Matrix notation denoted by [ ] has been used to simplify the above expression. For e.g., 
[7*^°)] denotes the 3-metric with components 7^°^ and [h] denotes hab, etc. In addition 
[71(0) (x)] is a matrix whose components are just the functional derivative of Q with 
respect to -f^^{x), 

sg SQ 



{19d) 
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where the last equahty follows from eq.(18&). [7f(o)] is its traceless counterpart where 
7r(o) denotes the trace: 

^ „,(0) ab —ab ^ab ^ ^ab ^ —(0) ^,(0) ^,(0) cd /or\\ 

no) = lab 7r(o) , 7r(o) = 71(0) - 3 7(0) 7r(o) , T^ab = lac Ibd 7r(o) ■ (20) 

In eq.(19c), the exponential of a matrix [A\ is defined in terms of a Taylor series 
expansion, 

exp[A] = [/] + [A] + i[A][A] + ... (21) 

At this point it is important to remember that \x{'^)Tlab{x)\ and [l^ab 

spatially dependent. Eq.(19a-d) represents a general solution to the strongly-coupled 

system consisting of gravity and dust. 



2.4.- Relation to Kasner Metric 

Using semiclassical methods, Francisco and Pilati (1985) have shown how to derive the 
Kasner metric describing pure gravity (see also Salopek and Stewart 1993). Here, it 
will be shown that the case of dust interacting with gravity admits a solution where the 
metric behaves like a Kasner metric at two separate epochs. 
(7/7(0))"'^^^ is a quadratic polynomial in x, 

with roots (x(i),X(2)), 




X{i) = \ 7(0)' m + Vo [7(o)7itoVi°6 



1/2 



(226) 



X(2) 



-1/2 
7(0) ^(0) 



■^-l=a6 :=(0) 
7(0)^(0)^a6 



1/2 



(22c) 



which are spatially dependent: X(i) = X(i)(^) ^'^'^ X(2) = X{i){^)- The average of the 
two roots. 



X(l) + X(2) 



O 1/2 

2 7(0) 7r(o) 



^—ab — (0) 



2 (^rfo) -67'(0)"a6 

gives the time when the extrema of (7/7(0))^^^ is reached, z becomes 



(23) 



3 - \ 1 _ ^ 



(24) 
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Note that the equation for [h] may be put in the following form 



[h] = [h\l^^] exp < 



2^ bm] b(o)] [7(of ] 



,(0)\V2 



[h 



l/2n 
(0)J 



(25) 



where the argument of the exponential is explicitly a symmetric matrix. Following 
Salopek (1992), we diagonahze the argument using an orthogonal matrix [O], 



where 



(26a) 



[Oy [O] = [O] [Oy = [/] (identity matrix) , 



(266) 



[D] = Diag[di,d2,4] , di + d2 + d3^0, di + d'^ + di^l. (26c) 
[h] then admits the following simple form: 

iT 



and [7] 



with 



[h] = [hl^] [O] 
Diag 



1 - 

X{2) 

1- ^ 

X(l) 









X 


X(2) 




X(2) 


X 1 






X 








X(l) 



[o] [hll^] 



[7] = [7(0^)1 [O] 



Diag 



1- — 

fi- — 
fi- — 

[O] [7(0^;] , 



(1) 



2p\ 



-) 



2pW 



2'p\ 



(3) 



1 - 



X(2)) 



(27) 



(28) 



(29a) 
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Given a spatial point for times x close to X{2){x), the diagonal part of the 3- metric 
in eq.(28) evolves like a Kasner universe, with Kasner exponents, p^\p^\p^^ ■ 

p^^^+pr+pr=i, {p^^^y+{p^^^)\{p'^^y=i. (29^) 

For times close to X{i){^)j the 3-metric also evolves like a Kasner universe, but with 
Kasner exponents, p^\pi'\p^^- (See Belinski et al 1970 and Salopek and Stewart 
1993.) 



2.5. Constant Functionals 



Given a generating functional S describing a HJ flow, C['~fab{x) , x{^)] is a Constant 
Functional if for all choices of the lapse N and shift iV\ C vanishes: 



I' 



SC . , , SC . / N 

x{x) + T — r-^7ab{x) 



6C 



6x{x 



(^N + N%) + 



(30) 



2 7V7-V2 [2jac{xhbd{^) - lab{xhcd{x)] 

where the last equahty follows from the HJ flow equations, 
(x - N%) /TV = 1 , 

{iab - Na\b - Nbla) /N = 2^'^^^ {2-faclbd - lablcd) 



+ Na\b + 



5C 



5S 



hcd{x) 



Since N , Ni are arbitrary, the constant functional C obeys the following: 

8C . , , , , , , , 5S SC 

+ 2 [2-/ac{x)%d{x) - ^ab{x)^cd{x 



5x{x) 



-2 7: 



0, 



ik 



6C 



One may show that if 



5C 5C 

T lkl,i + ^ X,i = . 

oiki ox 



Slab{x) j ' 

(31a) 
(31 &) 

(32 a) 
(326) 

(33) 



C^C[7r(x)] 

is an arbitrary gauge-invariant functional of the initial metric '^f^{x) defined by the 
minimization prescription eq.(186), then C is indeed a constant functional. Two simple 
examples are. 



j d^x R{Q) , etc. 



(34a) 
(34&) 
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Eq.(p3D is an elegant result which has wide ranging implications. For example, a 
complete list of constant functionals could be used to define a particular universe within 
the ensemble of universes whose evolution is described by the generating functional S. 
In fact, a space-time transformation of that particular universe would not change the 
numerical values of the complete list of constant functionals. 



3. Semiclassical Evolution: Dust and Gravity 

Solving the HJ equations, the energy constraint and the momentum constraint, brings 
us one step closer to understanding the quantum theory of the gravitational field 
because the generating functional S, eq.(18a), may be interpreted as the phase of the 
wavefunctional in the semiclassical approximation, 

^[7a6(x),x(a:)]~e^^[^-(^)'^(^)]/\ (35) 

where Planck's constant h is assumed to be tiny. Explicit solutions for S = 
<S['yab{x),x{x)], are now discussed. 

Semiclassical evolution is found by inverting the classical evolution eqs.(19a-d) to 
express the initial metric, 7^°'* = 'J^b {'Jed, x) j iii terms of the original fields, and then 
substituting to obtain S in eq.(18a). Typically, the inversion is very difficult because one 
must solve nonlinear, partial differential equations. Some special cases are given below 
and in Section 4 where the inversions may be performed explicitly without resorting 
to an approximation method. In the section denoted General Case below, one obtains 
semiclassical solutions by inverting the classical evolution equations through a Taylor 
series in x{^)- 



3.1. Elementary Example 

In the special case where the initial functional Q is proportional to the volume of the 
initial 3-geometry, 

Q = C J d^x 'j^Q^ , (where C is a homogeneous constant) (36) 

it is easy to perform the necessary inversion to compute S. Since tt^q-) vanishes, one finds 
that z = and that 

haJ = Kb , (37a) 

(1 - px) 

The generating functional S is then given by eq.(18a) 

SMx),x{x)] = -\l d'x^'/' ^ (3^ 
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This solution had be given previously by Salopek and Stewart (1992), and it is of the 
form given in eq.(|l]). It is mentioned here to show that the new formalism of this paper 
encompasses the earlier work. 

If C is real, one may define the classical Hubble parameter, H, through 

H{x) = -^^ab^, (39) 



ab 



giving 



2 1 
H{x) = - -, 

3 x(2;) - 



4 

3C 



(40) 



Hence, if C < and x{x) > 0, the Universe is expanding (locally) at each spatial point, 
whereas for C > and x(^) ^ 0? it is contracting locally until it reaches a singularity 
at x{^) = 4/(3C). The appearance of the singularity is problematic for the classical 
theory. C = yields the trivial solution S = 0. 



3.1.1. Evolution Beyond a Singularity Complex values of the parameter C are of 
interest because they help to explain how a universe may evolve beyond a singularity. 
Let us write 



■— = a — ib , where a, b are real , 
and then split S into real and imaginary parts: 
S = Sr + iSi , 



d^X'j 



1/2 



(x(x) - a) +52 



(41) 

(42 a) 
(426) 



46 

Sj=- I d'x-i 



,1/2 



(42c) 



(x(x) - a) +52 

The wavefunctional then becomes 

^ = exp {-Si/h) exp {iSR/h) . (43) 

As a result, one can define a classical Hubble parameter using the real part of the 
generating functional, Sr, in eq.( 

2 \x{x) 



H{x) 



{x{x) - a) +62 



(44) 



(Here I am assuming that H{x) or some closely related function can be represented 
by some Hermitean operator. In ordinary quantum mechanics, the expectation value 
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< p > of some Hermitean operator p is guaranteed to be real. Hence for a wavefunction 
which is tightly peaked about a classical trajectory, the expectation value of p may be 
computed approximately using the real part of the phase.) 

For [x{x) — a] large and negative in eg. (PD , a universe is contracting (locally), 
whereas for [x{x) — a] large and positive, the same universe is expanding (locally). 
Apparently, that universe was able to pass through the singularity at xi.^) = 0.1 and 
bounce from a contracting phase to an expanding phase. If 6 > 0, universes with large 
volumes are exponentially suppressed: 



|^'P = exp<(-^ I £x'j^^'^ 



1 



62 



(45) 



The above analysis was valid for arbitrary x{^)- These qualitative results concerning 
the bounce from a singularity are in agreement with a quantum formulation given by 
Salopek (1992b) where he assumed that x{^) was homogeneous. 



3.2. Intermediate Example 

The case of where the initial functional Q is a function of the initial volume is partially 
tractable: 



Q = f 



V, 



(0) 



with Vf 



(0) 



After minimizing with respect to 'j'^hix), the classical evolution equations give 



7(0) 



(46) 



Z = 0, hab{x 



/^i?(x) 



(47a) 
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1/2 



X 



1/2 
^(0) 



(476) 



where /'[V{o)] denotes the derivative of the function /[V(o)] with respect to the variable 
V(o). Hence, only the conformal factor of the 3- metric evolves. Bringing the conformal 
factor to the other side, and integrating over all x, one finds that, given a 3-metric 'jabi^) 
and field configuration x{^)^ ^(0) is defined implicitly through 



(0) 



£x^'/' - 



1 - 



2 ■ 



(48 a) 



The resulting generating functional for arbitrary x(x) is 
4 /• I/O 1 



5 = --/rf^xy/^ 



X - 



3/'[V(o)] 



(486) 
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By invoking a Legendre transformation, 



df 

.^/-^V^(o), (49a) 
^ (496) 



one can simplify eqs.(48a-b) even further, 

36 j 

fd'^^V^- 1 ^. (50&) 

One verifies that the second eq.(506) is a consequence of minimizing the first eq.(50a) 
with respect to the single independent variable h. Finally, using one more set of 
substitutions, 

4 

« = ^ , j(«) = 9{b) (51) 
one can express the generating functional in its simplest form, 

S[^ab{x),x{x)] d'x^'/' —l-—+j(a) , (52a) 

3 J ' [x{x) - da ^ ' 

Once again, the second eq.(52&) is obtained from the first by minimizing S with respect 
to a where j{a) is an arbitrary function of a. 

One could have derived this result more simply by starting with the solution, 

r[T*xWH = -i/Ay/^^^ (63) 

of the energy constraint and the momentum constraint, and by then constructing another 
solution S by taking a superposition over the homogeneous parameter a, 

Sbab{x),x{x)] = r[-fab,x{x)\a] +j{a) , (54a) 

which means that one chooses a according to the minimization prescription: 

da ^ da ^ ^ 

This solution was actually suggested by Salopek and Bond (1990) who gave a much 
different derivation. Here the aim was to show that the Green's function method 
encompasses previous results. In addition, the Green's function method has the 
advantage that it yields the classical evolution equations eq.(47a-b). 
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3.2.1. Explicit Computation for Intermediate Example If the arbitrary function j{a) 
is assumed to be linear, 

j(a) = a , (55a) 

one may readily compute the generating functional eq. (52a) , if x{x) = x is homogeneous, 
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.1/2 



(exact) , 



since one may invert eq.(526) to give. 



1/2 



(exact) . 



(55&) 



(55 c) 



The sign before the square root in eq.(55c) is arbitrary, and I have chosen it to be 
positive. (Choosing a negative sign yields another solution of the HJ equation.) 

One may obtain an approximate form which is valid for mildly inhomogeneous fields 
x{x), by assuming that a is large, and then expanding the implicit equation (52&) in 
powers of Ax(x) / a where 



^X{x) =x{x) -X, X = ^ J d^x'j^/^xix), V = j d^x^^'^ 
The first few terms of the implicit equation (52 &) yield. 



UV 



i + ^<(x-x)^ > 



+ . 



The generating functional then becomes 
with 



i6y 3 , ,2 



<(X-X)' >-\^ Jd'xj'/' {x{x)-xf ■ 
If x{x) = X is homogeneous, then one recovers the earlier result, eq.(55&). 



(56) 



(57) 



(58a) 



(586) 



3.3. General Case 

It is straightforward to write down the generating functional for S in terms of x(x) and 
the initial 3-metric: 



-I 



d'x 7(-o;/^ X{x) [2lf$,, Tfl? - \ vrfo) } + Qbi?] , 



5Q 



(59a) 



(596) 
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Using the roots (x(i),X(2)) of (7/7(0))^''^ defined in eqs.(22&-c), one may rewrite the 
generating functional in the elegant form: 

'5 = -^/ci^-7f/;-^+Q[7l?]. (60) 

J ^ ' X(l) X(2) 

For some applications, one would like to express this as a function of the original 
3-metric 7a6(a;). In general, the necessary inversion of eqs.(19a-19c/) is very difficult. 
However, in the limit as xi-^^) ~^ 0, this inversion is straightforward since 

^ab{x) 72V) ^ X{x) . (61) 
Hence from eq.(59a), one concludes that 

S ^ Q as x{x) ^0. (62) 
which justifies the claim made earlier in eq.(18(i). 

3.4- Advanced Example 

In order to be concrete, I will illustrate the main features of the General Case by 
assuming that 

jd^x 7^/)' [C + ER^o)] , (63a) 

in which case 

- ^ + ^^(0)) and 7(0^^(0) = -ER""' . (636) 



7(0) ^(0) 



3.4- 1- Early Time Behavior Expanding the classical solutions (19a-c) in a Taylor 
series in x(^)) finds 

o))'''' = l-^(3C + i?(o))x + ---, (64a) 



^ = 2£;V<)^i?X + ---, (646) 

[h]^[h^'^]{l-AE{R^o]b'-'^]x + ...) , (64c) 
where [R{o)] denotes the traceless Ricci tensor with contravariant indices: 

[R^r ^ Rt) - ll^) Rio) . (65) 

Only first order is shown here although one could expand to arbitrary order. Using an 
iterative process, one can invert these equations to obtain as a function of x and 7, 



TaA^) = lab{x) + X{x) {C^abix) + E [ARab{x) - R{x) -fab{x)]} + ... (66) 
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Substituting into eq.(59a), S becomes 

5 = y d^x'y^^^ {C + ER) + 



i - 2R'^Rai)j + ]^CER + -C^ 
This agrees with a similar calculation by Salopek (1997). 



3 
-( 

4 



+ 



(67) 



3.4-2. Late Time Behavior In order to determine the late time behavior of S , one 
introduces a factor of 7^/^ in eq.(59a): 



5 = jd'x^'/'x 





^^(0) ^ab 


3 ^(0) 




(l-ix7(ojNo))' 


2 A 7(0) ^(0) ^ab 



+ QbT], (68) 



In the limit that x is large and positive, S assumes the very simple form, 

S[^^{x),x{x)] = -\j d'x^'/' -ly . (69) 

However, it may not always be possible to take the limit to large x because one may 
first reach a singularity at a finite value of x- For this case, the form of the generating 
functional in terms of the original variables has not as yet been determined. 



4. Strongly-Coupled Solutions: Two Dust Fields Interacting with Gravity 

A non-trivial example which illustrates the main features of the Green's function solution 
method is provided by two dust fields, xi X2, interacting with gravity. The 
constraint equations for the long-wavelength system are: 

n^^\x)/K = r''' {'^laclU - 7a.7c.) + |^ + |^ = , (70a) 



'Hf\x)^-2\y,, 



hoh S-fcd Sxi Sx2 



5S \ 5S SS 6S 

7 + 7 lkl,i + T— Xl,i + y- X2,i = . 

^IkjJ O-fkl OXl 0X2 



After defining the average field, x{x)-i and the static field, y{x), through 
^ iX2 + Xl) ^ (X2 - Xi) 

X ^ , y 2 ' 

the constraints become 

n^'\x)/K = (27„,7,^ - 7„,7,^) + ^ = , 



Hf)(x) = -2(7 



S^ab S'^cd Sx 



6S \ SS 6S SS 



(70&) 

(71) 
(72 a) 
(726) 
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The static field y{x) appears in the momentum constraint eq.{72b) but it is absent in 
the energy constraint eq.(72a) which is identical to that of a single dust field, eg. (p!OD . 
For all intents, one can describe the two-dust system with gravity, which consists of the 
fields {xi^X2, lab), as a system containing a single dust field interacting with gravity, x 
and Jab, and containing a static field y{x) which appears in the generating functional S 
but which does not evolve. Hence given the initial functional, 

Q^Qb!S{x),y{x)] (73) 

which is invariant upon reparametrization of the spatial coordinates, 

the generating functional at later times is given by 

Sbab,x{x),y{x)] = g[-fab{x),x{x)\labi^)^X{o){x) = 0] + Q^l? (a^) , (x)] (75a) 

where Q is just the Green's function for dust with gravity defined in eq.(lla) with 
X(o){x) = 0, and 'y'^J{x) chosen to minimize Q + Q: 

^Tab hab 

The minimization prescription leads to the classical evolution equations (19a-d) which 
carry over directly to the present situation. 



Exact Solution for Two Dust Fields with Gravity 
Explicit solutions arise when the initial generating functional is of the form 



(TX 



1/2 



A + (t^'o) y,ay,b 



1/2 



(76) 



where A is a constant. Using the Green's function method described above, I will show 
in section 4.2 that the resulting generating functional S at arbitrary x{^) is 



S[iab{x),x{x),y{x)] 



d^x-y 



1/2 



+ 



d^x'y 



1/2 



1/2 



(77) 



This solution is a special case considered by Salopek, Stewart and Parry (1993). One 
may also demonstrate the validity of this solution by direct substitution into the HJ 
eq.(72a). 
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4-2. Derivation of Exact Solution 

The explicit solution (77) for two dust fields interacting with gravity will be justified 
using the Green's function method of Section 2. 

The momentum conjugate to 7^°^ is given by eq.{19d): 

1 



where 



s = 7fo) y,ay,b and |/'" = 7^0' y,b ■ 



ah 



(78a) 



(786) 



This example is qualitatively different from the exact solution given in Section 3 in that 

(79a) 



the traceless part of tt^q) is non- vanishing, 



-1/2 ab 
7(0) 7r(o) 



2sV2 



while the trace is given by 

-1/2 

7(0) ^(0) = ^ + s / ■ 
The square root of the determinant 7 evolves as a quadratic function of x, 



(7/7(0)) 



1/2 



1 - 



X(i). 



1 - 



X 

X(2), 



with roots at 



X(i) 



(796) 



(80a) 



(806) 



X{2) 



3A 



(80c) 



Please note that X{i) 7^ Xu X{i) denotes the first root of (7/7(0))"^^^ whereas Xi denotes 
the first dust field. The 3-metric evolves according to 

4/3 



7ab = 1 



X 

X{2), 



B. 



ac 7(0) 



(81a) 



where the symmetric matrix Bac has components 



^ac lac ^ 



_X 

X(2 ) 

1- ^ 

X(2) 



?/,a y,c 



(816) 



The inverse 3-metric has components 

-4/3 



7 



ab 



X 

X{2), 



/^ac „,(0) /^db 
7cd • 



(82 a) 
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where C"''^ is a symmetric matrix: 



^ac ^ ^ac _^ 1 V^W ^;a^;c 



1 - 

X(i)' 



(826) 



One concludes that 



\ 2/3 

1 - ^ 

l/|a = 7'^" 2/,a y,t = \ ""'''^2 7(0) y,a V • 



1 



X(i) 



(83) 



One may readily invert this relation to find s, eq.(78&), as function of the original fields: 

(84) 



.1/2 





(1- 








(1 - ¥x) 




-1/2- 



Using eg. (|60D , one can write S in terms of x(x) and the initial 3-metric, 



S= Id'x^li^ 



15) 



Eliminating the initial variables in favor of the original ones, one recovers the claimed 
result eq.(77). 

5. Strongly-Coupled Solutions: Scalar Field Interacting with Gravity 

5.1. Massless Scalar Field, Gravity and Cosmological Constant 

The strongly-coupled system of gravity interacting with a scalar field is also tractable. 
Most of the techniques that applied to dust and gravity can also be applied here although 
the details are quite different. The main results will be stated without much elaboration. 

In the strongly-coupled limit, the Hamiltonian density for gravity and a scalar field 
with cosmological constant is: 

2 



A complete solution ("Green's function solution") for the generating functional 
describing a massless scalar field with cosmological constant term interacting with 
gravity is: 



Qbab{x),(t){x)\ 7af, (a;),0(o)(a;)] = — J d X 

7 + 7(0) - 27^/^7(0^)^ cosh ( \/f \/^^ + ~ 



1/2 



(87) 
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where z was defined in eq.(ll&). The sign of the Green's function is arbitrary and a 
minus sign was chosen in eq.(87). Hereafter, I will assume that 0(o) = for similar 
reasons that were given for a dust field in Sections 2.1 and 2.2. 
In order to satisfy the momentum constraint, 

(bS \ bS 

one constructs a solution S using the Superposition Principle: 

S[jab{x), <P{x)] = Gbabix), m\ labi^), <P(o){x) = 0] + Q^i?] , (89a) 
where has been chosen to minimize Q + Q, 

SG sa 



Q is an arbitrary gauge-invariant functional of 7^^'' . Please note that S coincides with 
Q when = 0: 

S[yab{x),(l){x) = 0] = Q[^ab{x)]. (90) 

5.2. Classical Evolution 
Classical evolution is given by: 



z — (f) 



3 '(0)^(0) ^ ^0 



0)^(0) ^ab 



1, 



(916) 



[h] = exp 



2^ [W(o)] [7^°)] 

^(0) ^ab 



(91c) 



A and 9 are simply abbreviations for the following expressions: 

^ -1/2 

A 



/(()) 



i ((}) 



(92a) 



o7(o) ^('o) - H / \/ o7(o) ^('o) - 27(0)^(0) <b ' Vo , 



1^2 



-lr=ab ^(0) 



(926) 



6. Semiclassical Evolution: Scalar Field, Gravity and Cosmological 
Constant 



23 



It is straightforward to write S in terms of and the initial 3-metric 

sinh^ 



VA^ - 1 cosh6' + y4sinh6' 



(93) 



with A and 6 defined in eqs.(92a,b), but it is very difficult to find an explicit expression 
in terms of the original variables. I will be content to illustrate the 'early' and 'late' 
time behavior. 



6.1. Behavior for small (j){x) 

To this aim, note that for small the original 3-metric evolves according to 

V2 7(-o;/' 



(0) , 

7afe = Tab + 



This may be inverted by iteration to give Tab (x) as a function of and 7afe(a;): 



(0), 



9„(0) (0)' 
^^ab - ^(0) lab 



+ 



(94) 



(0) _ 
lab ~ lab 



7 



-1/2 



h-1 7r2 - 27-i7f'^*' Wab - Vo 



[2nab - IT -fab] + ... . 



(95) 



Substituting into eg. (p3D , one finds to first order in (f){x) that 
S['yab{x),(j){x)] = QlTbix)] + 

V2 I d^xj'/^ ^^7-1 7r2 - 27-i7f«^ 7rab-Vo + ... 
(small 0(x) behavior) 



(96a) 



with 



TT 



ab 



6Q 



ab 



_ _afe r^^A -^oh „ab „ -,ab 

, TT = 7af, TT , and vr - ^ 



-TT 7 



(966) 



which is in agreement with Salopek (1997). 
6.2. Behavior for large '^'^/'^{x) 

In order to determine the behavior at large 7^*^^, one rewrites the generating functional 
S as 



/ \ V2 

l + ^-2m) cosh^ 
7 \ 7 / 



1/2 



+ Q[7l?]- (97) 
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Provided — 1 > 0, 7 becomes very large as ^ — > Ocru where 9crit is given by 
(^crit = -tanh . 



(98) 



In this hmit, 7(o)/7 0, and S is proportional to the volume of any given 3- geometry, 



S[jabix), Hx)] = ^ 1 d'xj'/^ , (large j^/^ix) behavior) . (99) 
For the sake of brevity, other cases will not be considered here. 



7. Massless Scalar Field and Gravity 

A complete solution for the generating functional describing gravity and a massless 
scalar field without cosmological constant is: 



GMx), 0(a;)| 7l?(a;), 0(o)(a;)] = - d'x 7^/' 7(o^)' exp-J ^z^ + {^-^^,^y ,(100) 



where z is still defined by eq.(ll&). I will assume again that 0(o) = 0. 

7.1. Classical Evolution 
Classical evolution is given by 



V7(o)y 



1/4 



-1/2 



-7(0) 7r(o)exp 



7^(0) 



(0) 



(0) "a& 



(101a) 



(1016) 



[/.] = [/.(°)] 



exp 



2^ [%)] [7'°>] 



(0)\V2 



In particular for 

Q = / d^^x 7;/; 
the classical evolution is given by. 



C + ER^o) 



'-) 

.7(0)7 



1/4 



2(C + |i?(o)) exp 



(101c) 



fl02a) 



1 - 



SE^^^b ^(0) 
3 -'l'(0)-'l'afe 

(C + f i?(o)) 



l/2\ 



,(102&) 



25 



3 C 



(0) 



1/2 



- 1 



(102c) 



[/i] = exp 



[^(0)] b^'n 



^J(C+ f i?(o) ) - ^i?"o)i?l°J 



(102d) 



Semiclassical evolution is very similar to the case with cosmological constant, and 
it will not be discussed further. 



8. Conclusions 

In the semiclassical approximation, it was shown how to solve the system of strongly- 
coupled gravity and matter using a Green's function method. There are two steps in 
implementing this program: 

1. Computing the Green's function. Explicit Green's function solutions of the energy 
constraint in a Hamilton-Jacobi context were given for systems describing Einstein 
gravity interacting with either a dust field or a scalar field. 

2. Constructing the General Solution through the Superposition Principle. One can 
construct a general semiclassical solution to the momentum constraint as well as the 
energy constraint by a superposition over the parameter fields of the Green's function. 

The number of parameter fields appearing in the Green's function is a tricky issue. 
Apparently, this number should be one less than the number of fields originally in the 
problem. For example, in the case of a dust field x{^) interacting with the 3-metric 
7a6(,T), there arc 7 degrees of freedom per spatial point x: one for the dust field and 6 
for the symmetric 3x3 matrix 'fab- Hence, there should be at most 6 parameter fields in 
the Green's function. If there are more, one should set the additional parameter fields 
to zero. Exact solutions demonstrate the validity of this approach. 

The problem of strongly-coupled gravity interacting with matter is not exceedingly 
complicated, although the actual details may be intricate. Many of the exact solutions 
presented in this paper had been derived earlier using other techniques. However, the 
Green's function method is very general and in some sense all solutions to strongly- 
coupled gravity and matter may be derived using it. The general methodology may 
perhaps be useful in analyzing a wide range of gravitational phenomena. 
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